An optimized chiral nucleon-nucleon interaction at next-to-next-to-leading order 
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We optimize the nucleon-nucleon interaction from chiral effective field theory at next-to-next- 
to-leading order. The resulting new chiral force NNLO op t yields x 2 ~ 1 P er degree of freedom for 
laboratory energies below approximately 125 MeV. In the A — 3, 4 nucleon systems, the contributions 
of three-nucleon forces are smaller than for previous parametrizations of chiral interactions. We use 
NNL0 O pt to study properties of key nuclei and neutron matter, and we demonstrate that many 
aspects of nuclear structure can be understood in terms of this nucleon-nucleon interaction, without 
explicitly invoking three-nucleon forces. 



PACS numbers: 21.30.-x, 21.10.-k, 21.45.-v, 21.60.-n 

Introduction - Interactions from chiral effective field 
theory (EFT) employ symmetries and the pattern of 
spontaneous symmetry breaking of quantum chromody- 
namics [TJ [2]. In this approach, the exchange of pi- 
ons within chiral perturbation theory yields the long- 
ranged contributions of the nuclear interaction, while 
short-ranged components are included as contact terms. 
The interaction is parametrized in terms of low-energy 
constants (LECs) that are determined by fit to exper- 
imental data. The interactions from chiral EFT ex- 
hibit a power counting in the ratio Q/A, with Q being 
the low-momentum scale being probed and A the cut- 
off, which is of the order of 1 GeV. At next-to-next-to- 
leading order (NNLO), three-nucleon forces (3NFs) enter, 
while four-nucleon forces (4NFs) enter at next-to-next- 
to-next-to-leading order (N 3 LO). For laboratory energies 
below 125 MeV, the nucleon-nucleon (NN) force exhibits 
a quality of fit with x 2 ~ 10/datum at NNLO while a 
high-precision potential N 3 LOem 5 with a x 2 ~ 1/datum 
up to 290 MeV, was obtained by Entem and Machleidt 

Hi- 

The 3NFs at NNLO that accompany the current N 3 LO 
NN potentials play a pivotal role in nuclear structure 
calculations [5] . They determine the ground-state spin of 
10 B [6], correctly set the drip line in oxygen isotopes [HIE], 
and make 48 Ca a doubly magic nucleus (9J, [10] . While it 
might seem surprising that smaller corrections at NNLO 
are so decisive for basic nuclear structure properties, the 
3NF contains spin-orbit and tensor contributions that 
clearly are important for the currently employed chiral 
interactions. The contributions of 3NFs at N 3 LO have 
also been worked out [H] [12], and there are on-going 
efforts to compute even higher orders [13] . 



While the quest for higher orders is important, this ap- 
proach will result in higher accuracy only if the optimiza- 
tion at lower orders was carried out accurately. Thus, it 
is important and timely to revisit the optimization ques- 
tion. We note in particular that the fits of the currently 
employed chiral interactions [3] H] [14] date back about 
a decade and that there has been a considerable recent 
progress in developing tools for the derivative-free non- 
linear least-squares optimization [15] . Furthermore, the 
quantification of theoretical uncertainties is a long-term 
objective of nuclear structure theory, and this requires 
a covariance analysis of the interaction parameters with 
respect to the experimental uncertainties of the nucleon- 
nucleon elastic scattering observables; see, for example, 
Refs. [IS] [TB]. This letter takes the first step toward 
this goal. We present a state-of-the-art optimization of 
the NN chiral EFT interaction at NNLO. This yields 
a much-improved x 2 arL d a high-precision NN potential 
NNLO opt . The 3NF at NNLO is adjusted to the bind- 
ing energies in A — 3, 4 nuclei. We present computations 
of three-nucleon and four-nucleon bound states, and we 
employ NNLO op t to ground states and excited states in 
10 B, masses and excited states of oxygen and calcium 
isotopes, and neutron matter. 

Optimizing the NN interaction at NNLO - For the op- 
timization of the chiral AW interaction we use the Prac- 
tical Optimization Using No Derivatives (for Squares) al- 
gorithm, POUNDerS [T5], as implemented in [17]. This 
derivative-free algorithm employs a quadratic model and 
is particularly useful for computationally expensive ob- 
jective functions. We optimize the three pion-nucleon 
(irN) couplings (01,03,04), and 11 partial wave contact 
parameters C and C, while we keep the axial-vector cou- 
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pling constant gA, the pion-decay constant /„■, and all 
masses fixed. In the optimization, we minimize the ob- 
jective function 
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where <5 NNLO are NNLO phase shifts, S^™ 93 are ex- 
perimental phase shifts from the Nijmegen multi-energy 
partial- wave analysis [18], x denotes the parameters of 
the chiral interaction, and w q are weighting factors. Note 
that Eq. ^ is not the x 2 with respect to experimen- 
tal data. The actual \ 2 is calculated following the 
POUNDerS optimization. The phase shifts <5 NNLO are 
computed from i?-matrix inversion, and in the proton- 
proton (pp) channels we include the Coulomb interac- 
tion [113 [5D]. The contact terms are optimized to re- 
produce the Nijmegen phase shifts for each correspond- 
ing partial wave, while keeping the Cj's fixed. For the 
contacts, the weight w q scales with the third power of 
the relative momentum q, while for the Cj's, we em- 
ploy the uncertainties quoted in the Nijmegen analy- 
sis [TS]- This approach can be justified by a physi- 
cal argument: for the peripheral waves the higher en- 
ergies still represent longer-range physics, and the need 
for a pedantic agreement with lower energy phase shifts 
can be weakened. The irN couplings Ci,c 3 , and c 4 
were simultaneously optimized to the peripheral partial- 
waves 1 D 2 , 3 D 2 , 3 F 2 ,E 2 , 3 F 3 , 1 G 4l and 3 F 4 . Note that 
the NNLO contact terms do not contribute to orbital 
angular momenta L > 2. We do not include other pe- 
ripheral waves from the Nijmegen study since they carry 
extremely small uncertainties, which lead to a very noisy 
objective function. 

Table [I] summarizes the optimization results. Our val- 
ues should be compared with the irN couplings as deter- 
mined from irN scattering data, where c% = — 0.81±0.15, 
c 3 = -4.69 ± 1.34, and c 4 = +3.40 ± 0.04 have been ob- 
tained 21j. Thus, POUNDerS yields values for c 4 and C3 
that agree well with the empirical determination from nN 
scattering. The c 4 value, however, deviates significantly 
from its empirical value. The same trend was found in 
the construction of the N 3 LO [I] NN interaction. A 
detailed statistical sensitivity analysis of the LECs with 
uncertainty quantification will be presented in Ref. |22j . 

Table |n] shows the x 2 /datum for NNLO op t at various 
laboratory energy bins. The quality of the fit is particu- 
larly good for energies below 125 MeV. For comparison, 
the np NNLO interaction of Ref. [3] yields % 2 /datum of 
12-27 in the range A = 600/700 - 450/500 MeV at ener- 
gies up to 290 MeV. 

Around energies of 144 MeV there exist two data sets 
of pp differential cross sections with a very high precision 
(0.5% error) [Mj (47 data points). The total number 
of pp data in the energy interval 125-183 MeV is 343. 



TABLE I. Parameters of NNLO opt at A = 500 MeV and 
Asfr = 700 MeV [3]: c, (in GeV" 1 ), C (in 10 4 GeV" 2 ), and 
C (in 10 4 GeV" 4 ). 



LEC 



Value 



I LEC Value 



LEC 



Value 



Cl 

Cl s 
Ci Pl 

C3 Sl -3D 1 



-0.91863953 
-0.15136604 
2.40402194 
0.41704554 
0.61814142 



C3 
C3 Sl 

C 3 p 2 



-3.88868749 
-0.15214109 
0.92838466 
1.26339076 
-0.67780851 



c 4 

/-inn 

^1 S 
£3 Si 



a 



Pi 



4.31032716 
-0.15176475 
-0.15843418 
-0.78265850 



TABLE II. x 2 /datum for NNLO opt at A = 500 MeV and 
Asfr = 700 MeV [3] with respect to the np and pp 1999 
databases |23| . The values without the high-precision data 
sets [24] are marked by asterisks. 



T lab (MeV) 


0-35 35-125 125-183 183-290 


0-290 


pp x 2 /datum 
np x 2 /datum 


1.11 1.56 • 
0.85 1.17 


1.87 6.09 


\ 


17.10 
14.03* 
2.95 



The unusual precision of those 47 data points distorts 
the x 2 /datum for this interval. For this reason, Table [n] 
also shows the results without the high-precision data. 

Two comments arc in order. First, the x 2 with respect 
to scattering observables is lower when the 1 P\ phase 
shifts arc weighted with the uncertainties from the Ni- 
jmegen analysis. The P-waves are accurately reproduced 
only when going to N 3 LO [4]. Second, the 3 Si — 3 D\ cou- 
pled channel is optimized with the additional constraint 
of reproducing the deuteron binding energy. The remain- 
ing deuteron observables, as well as the LSo scattering 
observables, are predictions and reproduce the experi- 



mental values well; see Table III 

Figure [T] shows some np phase shifts of NNLO opt and 
compares them with phase shifts from other potentials 
and partial wave analyses. Apart from the 3 P-waves, 
the phase shifts of NNLO op t closely agree with those ob- 
tained at N 3 LO. Note, however, that these deviations 
do not spoil the good x 2 a t laboratory energies below 
125 MeV. 

Three-nucleon forces also appear at NNLO, and two 
additional LECs (cd and ce) enter. These are deter- 
mined from calculations in the three-nucleon and four- 
nucleon systems. We find that the binding energies of 3 H, 
3 He, and 4 He do not uniquely determine cd and ce, and 
the parametric dependence of both LECs is very similar 
to those found in previous studies [32l [33] . Therefore, 
we choose c D — —0.2 guided by the triton half life [33 
and obtain ce = —0.36 from optimization to the binding 
energies. The resulting point charge radii of 4 He are also 



in good agreement with experiment; see Table IV 

Performance of NNLO opt for light- and medium-mass 
nuclei and neutron matter - In this paper, we apply 
NNLO op t to 10 B, isotopes of oxygen and calcium, and 
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TABLE III. Scattering lengths a and effective ranges r (both 
in fm). The superscripts N and C for the proton-proton 
observables refer to nuclear forces and Coulomb-plus-nuclear 
forces, respectively. Bd, td, Qd, and Pd denote the deuteron 
binding energy, radius, quadrupole moment, and D-state 
probability, respectively. Qd and m are calculated without 
meson-exchange currents and relativistic corrections. 



N 3 LOem NNLOopt Exp. Ref. 
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FIG. 1. (Color online) Computed np phase shifts of the opti- 
mized NNLO potential of this work (red), the NNLO potential 
of Ref. .3; (dashed, blue), and the N 3 LO potential [4] (green, 
dotted) compared with the Nijmegen phase shift analysis |18| 
(solid dots) and the VPI/GWU analysis SM99 [31] (open cir- 
cles) . 



neutron matter. The considered systems are particularly 
interesting because the current NN chiral interactions 
at N 3 LO completely fail to describe key aspects of their 
structure. 

To study the ground- and first excited state in 10 B, we 
carry out no-core shell model (configuration interaction) 
calculations [31] using the bare NNLO op t in model spaces 



TABLE IV. Ground-state energies (in MeV) and point proton 
radii (in fm) for 3 H, 3 He, and 4 He using the NNLO op t with 
and without the NNLO 3NF interaction for cd = —0.20 and 
c E = -0.36. 





£( 3 H) 


£( 3 He) 


.E( 4 He) 


r P ( 4 He) 


NNLO 


-8.249 


-7.501 


-27.759 


1.43(8) 


NNLO+NNN 


-8.469 


-7.722 


-28.417 


1.43(8) 


Experiment 


-8.482 


-7.717 


-28.296 


1.467(13) 



of up to iV max = 10 harmonic oscillator (HO) shells 
(10 hfl) above the unperturbed configuration. These 
model spaces are not large enough to provide fully con- 
verged results for the ground state and first excited state 
of 10 B. Still, the variational upper bounds for the ener- 
gies are -54.35 MeV for the 1+ state and -54.32 MeV for 
the 3 + state. The energies are very close, in contrast to 
N 3 LOem> which yields a level spacing of about 1.2 MeV 
between the J 71 ' = 1 + ground state and the J 71 " = 3 + 
excited state [BJ. 

Chiral NN interactions at N 3 LO fail to explain the 
neutron drip-line in oxygen isotopes, and 3NFs have been 
the key element for understanding the structure of nu- 
clei around 24 O [7JIE]- Figure [2] shows the experimental 
ground-state energies of oxygen isotopes and compares 
the results from coupled-cluster (CC) computations in 
the A triples approximation 35 -37 . Our CC calculations 
employ a Hartree-Fock basis (HF) built from iV max = 15 
HO shells at Ml = 20 MeV. Because of the "softness" of 
NNLOopt, this model space is sufficiently large to con- 
verge the ground states and excited states of the nu- 
clei considered. In addition, we performed shell-model 
(SM) calculations assuming the closed 16 core with an 
effective interaction derived from many-body perturba- 
tion theory to third order in the interaction and includ- 
ing folded diagrams [35]. For the SM calculations, the 
single-particle energies were taken from the experimen- 
tal 17 spectrum. In both CC and SM, NNLO opt results 
are close to experiment. In contrast, the N 3 LOem case 
requires 3NFs to provide reasonable description of mea- 
sured values. 

Now we consider the heavy isotopes of calcium. Here, 
48 Ca is doubly magic, 52 Ca exhibits a soft subshell clo- 
sure, and 54 Ca is predicted to have an even softer subshell 
closure 10J. A signature of shell closure is the location of 
the first 2 + state. We employed CC equation-of-motion 
methods within the singles and doubles approximation 
[571 [55] to compute the first 2 + state in the calcium iso- 
topes. Figure [3] shows that N 3 LOem fails to describe 
the location of the first 2+ state in 40,48,50,52,54,56 C£L 
In contrast, NNLO opt yields 48 Ca as a doubly magic 
nucleus and predicts subshell closures in 52,54 Ca. The 
NNL0 O pt overbinds the calcium isotopes by about 1 MeV 
per nucleon. In particular 40 ' 48 ' 52 Ca are overbound by 
1.03 MeV, 1.06 MeV, and 1.04 MeV per nucleon, respec- 
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FIG. 2. (Color online). The ground-state energies of oxygen 
isotopes obtained in CC with the NNLO op t and N 3 LOem in- 
teractions compared with experiment. The inset shows SM 
results. 



binding energies of oxygen and calcium isotopes. How- 
ever, excitation energies relative to the ground state con- 
verge somewhat better. Our results for the first 2 + state 
in 22 < 24 and in 48 Ca are 2.3(3) MeV, 3.5(5) MeV and 
4.8(7) MeV, respectively. We estimate the uncertainty 
by varying Ml in the interval 16-22 MeV. The results ob- 
tained by using NNLO op t NN interaction alone yields 
2.5 MeV, 5.0 MeV, and 4.5 MeV in 22 < 24 and 48 Ca, re- 
spectively. These preliminary results suggest that the 
3NFs may not dramatically change the results that were 
obtained with the NNLO op t NN interaction alone. 

It is instructive to compare the predictions of NNLO op t 
and N 3 LOem for the neutron matter equation of state 
at sub-saturation densities with the results of ab-initio 
calculations of Rcfs. [32]. Figure [4] shows that the per- 
formance of NNL0 O pt is on par with the EGM results of 
Ref. [32] , which take into account the effects of 3NFs and 
4NFs. The predictions of N 3 LOem deviate from other re- 
sults at higher densities. 
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FIG. 3. (Color online). The first 2^ state in selected calcium 
isotopes obtained in CC with the NNLO op t and N 3 LOem in- 
teractions compared with experiment. 
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tively. That is, the excess energy per nucleon is fairly con- 
stant; hence, NNLO op t reproduces binding energy differ- 
ences, such as neutron-separation energies and low-lying 
excited states, rather well. 

The complete description of nuclei at NNLO also 
requires 3NFs. We computed the first 2 + state in 
22 ' 24 and in 48 Ca with the 3NF compatible with the 
NNLOopt interaction. The matrix elements of the 3NF 
are expensive computationally, and we must at present 
limit their calculation to three-body energies up to 
e3max = 2n a + l a + 2n 6 + / 6 + 2n c + l c = 14. (Recall 
that we employ 15 major harmonic oscillator shells for 
the NN interaction.) We also used the normal ordered 
two-body approximation for the 3NF 10 41] with respect 
to a HF reference. With the restriction of e3 max = 14, 
we were not able to obtain fully converged results for the 



FIG. 4. (Color online). Energy per nucleon for neutron mat- 
ter for NNL0 O pt and N 3 LOem [I]- The calculations used 
the CC method with the inclusion of particle-particle ladders 
and a continuous single-particle spectrum. The shaded area 
(EGM) shows uncertainty bands for N 3 LO chiral effective field 
theory calculations of Ref. [42], including 3NFs. 

Conclusions - We constructed the new NN chiral EFT 
interaction NNLO op t at next-to-next-to-leading order us- 
ing the optimization tool POUNDerS in the phase-shift 
analysis. The optimization of the low-energy constants 
in the iViV-sector at NNLO yields a x 2 /datum of about 
one for laboratory scattering energies below 125 MeV. 
The NNL0 O pt NN interaction yields very good agree- 
ment with binding energies and radii for A = 3,4 nu- 
clei. Key aspects of nuclear structure, such as excita- 
tion spectra, the position of the neutron drip line in oxy- 
gen, shell-closures in calcium, and the neutron matter 
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equation of state at sub-saturation densities, are repro- 
duced by NNLOopt interaction alone, without resorting 
to 3NFs. We performed the initial calculation of the 
first 2+ states in 22 < 24 and 48 Ca with NNLO opt supple- 
mented by a 3NF and found effects of 3NFs to be small 
and good agreement with experimental excitation ener- 
gies. The precise role of 3NFs in medium-mass nuclei, the 
quantification of theoretical uncertainties, and optimiza- 
tions at higher-order chiral interactions will be addressed 
in forthcoming investigations. 
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